91 research outputs found
A mixed finite volume scheme for anisotropic diffusion problems on any grid
We present a new finite volume scheme for anisotropic heterogeneous diffusion
problems on unstructured irregular grids, which simultaneously gives an
approximation of the solution and of its gradient. In the case of simplicial
meshes, the approximate solution is shown to converge to the continuous ones as
the size of the mesh tends to 0, and an error estimate is given. In the general
case, we propose a slightly modified scheme for which we again prove the
convergence, and give an error estimate. An easy implementation method is then
proposed, and the efficiency of the scheme is shown on various types of grids
Quelques Résultats sur les Espaces de Sobolev
Nous établissons quelques résultats classiques sur les espaces de Sobolev, sur des ouverts à bord "faiblement" Lipschitzien. Ces ouvers sont ceux pouvant être localement ramené, par un homeomorphisme bi-lipschitizien, à un bord plat - ce qui est plus faible qu'être localement l'épigraphe d'une fonction lipschitzienne. Après avoir couvert les notions géometriques nécessaires sur ces ouverts (transport par des homéomorphismes bi-lipschitizeins, définition de l'intégrale sur le bord, etc.), nous donnons la définition d'un espace de Sobolev, et étudions les opérateur de prolongement, densité des fonctions régulières, injection de Sobolev et Rellich, trace etc
A Hybrid High-Order method for Leray-Lions elliptic equations on general meshes
In this work, we develop and analyze a Hybrid High-Order (HHO) method for
steady non-linear Leray-Lions problems. The proposed method has several assets,
including the support for arbitrary approximation orders and general polytopal
meshes. This is achieved by combining two key ingredients devised at the local
level: a gradient reconstruction and a high-order stabilization term that
generalizes the one originally introduced in the linear case. The convergence
analysis is carried out using a compactness technique. Extending this technique
to HHO methods has prompted us to develop a set of discrete functional analysis
tools whose interest goes beyond the specific problem and method addressed in
this work: (direct and) reverse Lebesgue and Sobolev embeddings for local
polynomial spaces, -stability and -approximation properties for
-projectors on such spaces, and Sobolev embeddings for hybrid polynomial
spaces. Numerical tests are presented to validate the theoretical results for
the original method and variants thereof
Convergence in of weak solutions to perturbed doubly degenerate parabolic equations
We study the behaviour of solutions to a class of nonlinear degenerate
parabolic problems when the data are perturbed. The class includes the Richards
equation, Stefan problem and the parabolic -Laplace equation. We show that,
up to a subsequence, weak solutions of the perturbed problem converge
uniformly-in-time to weak solutions of the original problem as the perturbed
data approach the original data. We do not assume uniqueness or additional
regularity of the solution. However, when uniqueness is known, our result
demonstrates that the weak solution is uniformly temporally stable to
perturbations of the data. Beginning with a proof of temporally-uniform,
spatially-weak convergence, we strengthen the latter by relating the unknown to
an underlying convex structure that emerges naturally from energy estimates on
the solution. The double degeneracy --- shown to be equivalent to a maximal
monotone operator framework --- is handled with techniques inspired by a
classical monotonicity argument and a simple variant of the compensated
compactness phenomenon.Comment: J. Differential Equations, 201
The gradient discretisation method for linear advection problems
We adapt the Gradient Discretisation Method (GDM), originally designed for
elliptic and parabolic partial differential equations, to the case of a linear
scalar hyperbolic equations. This enables the simultaneous design and
convergence analysis of various numerical schemes, corresponding to the methods
known to be GDMs, such as finite elements (conforming or non-conforming,
standard or mass-lumped), finite volumes on rectangular or simplicial grids,
and other recent methods developed for general polytopal meshes. The scheme is
of centred type, with added linear or non-linear numerical diffusion. We
complement the convergence analysis with numerical tests based on the
mass-lumped P1 conforming and non conforming finite element and on the hybrid
finite volume method
Unified convergence analysis of numerical schemes for a miscible displacement problem
This article performs a unified convergence analysis of a variety of
numerical methods for a model of the miscible displacement of one
incompressible fluid by another through a porous medium. The unified analysis
is enabled through the framework of the gradient discretisation method for
diffusion operators on generic grids. We use it to establish a novel
convergence result in of the approximate
concentration using minimal regularity assumptions on the solution to the
continuous problem. The convection term in the concentration equation is
discretised using a centred scheme. We present a variety of numerical tests
from the literature, as well as a novel analytical test case. The performance
of two schemes are compared on these tests; both are poor in the case of
variable viscosity, small diffusion and medium to small time steps. We show
that upstreaming is not a good option to recover stable and accurate solutions,
and we propose a correction to recover stable and accurate schemes for all time
steps and all ranges of diffusion
A unified analysis of elliptic problems with various boundary conditions and their approximation
We design an abstract setting for the approximation in Banach spaces of
operators acting in duality. A typical example are the gradient and divergence
operators in Lebesgue--Sobolev spaces on a bounded domain. We apply this
abstract setting to the numerical approximation of Leray-Lions type problems,
which include in particular linear diffusion. The main interest of the abstract
setting is to provide a unified convergence analysis that simultaneously covers
(i) all usual boundary conditions, (ii) several approximation methods. The
considered approximations can be conforming, or not (that is, the approximation
functions can belong to the energy space of the problem, or not), and include
classical as well as recent numerical schemes. Convergence results and error
estimates are given. We finally briefly show how the abstract setting can also
be applied to other models, including flows in fractured medium, elasticity
equations and diffusion equations on manifolds. A by-product of the analysis is
an apparently novel result on the equivalence between general Poincar{\'e}
inequalities and the surjectivity of the divergence operator in appropriate
spaces
A unified approach to Mimetic Finite Difference, Hybrid Finite Volume and Mixed Finite Volume methods
We investigate the connections between several recent methods for the
discretization of anisotropic heterogeneous diffusion operators on general
grids. We prove that the Mimetic Finite Difference scheme, the Hybrid Finite
Volume scheme and the Mixed Finite Volume scheme are in fact identical up to
some slight generalizations. As a consequence, some of the mathematical results
obtained for each of the method (such as convergence properties or error
estimates) may be extended to the unified common framework. We then focus on
the relationships between this unified method and nonconforming Finite Element
schemes or Mixed Finite Element schemes, obtaining as a by-product an explicit
lifting operator close to the ones used in some theoretical studies of the
Mimetic Finite Difference scheme. We also show that for isotropic operators, on
particular meshes such as triangular meshes with acute angles, the unified
method boils down to the well-known efficient two-point flux Finite Volume
scheme
Convergence analysis of a mixed finite volume scheme for an elliptic-parabolic system modeling miscible fluid flows in porous media
International audienceWe study a Finite Volume discretization of a strongly coupled elliptic-parabolic PDE system describing miscible displacement in a porous medium. We discretize each equation by a finite volume scheme whose properties are to handle a wide variety of unstructured grids (in any space dimension) and to give strong enough convergence to handle the nonlinear coupling of the equations. We prove the convergence of the scheme as the time and space steps go to . Finally, we provide numerical results which show that, from a practical point of view, the scheme behaves very well
- …